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Abstract 

Background: Evidence is increasingly accumulated about multiple roles for the p2-adrenoceptor gene in asthma. The results 
were inconsistent partly due to small sample sizes. To assess the association between p2-adrenoceptor gene 
polymorphisms and asthma risk, a meta-analysis was performed. 

Methods: We comprehensively searched the PubMed, EMBASE, BIOSIS Previews databases and extracted data from all 
eligible articles to estimate the association between (32-adrenoceptor gene polymorphisms and asthma risk. The pooled 
odds ratio (OR) with 95% confidence intervals (CIs) were calculated. 

/fesi//f5:Thirty-seven studies involving 6648 asthma patients and 1 5943 controls were included in the meta-analysis. Overall, 
significant associations were found in allelic genetic model (OR=1.06, 95% CI = 1 .01 — 1.1 2), recessive genetic model 
(OR = 1 .11 , 95% CI = 1.02-1.21) for Arg/Gly16. Stratified by ethnicity and age, significant associations were also found in 
Asian population in allelic genetic model, recessive genetic model and addictive model. For Gln/Glu27, no significant 
association was found when we combined all eligible studies. Age stratification showed significant associations in adults in 
allelic genetic model and recessive genetic model, but no significant association was found among Asians and Caucasians in 
ethnicity stratification. 

Conclusions: This meta-analysis implied that the p2-adrenoceptor Arg/Gly16 polymorphism was likely to contribute to 
asthma risk in Asian population. Gln/Glu27 polymorphism might be a contributor to asthma susceptibility for adults. 
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Introduction 

Asthma is a common chronic disorder of the airways, which is 
characterized by airway hyperresponsiveness, obstruction, and 
airway wall remodeling [1]. Based on data from World Health 
Survey, To T and colleagues found that the global prevalence of 
doctor-diagnosed asthma in adults was estimated at 4.3% [2]. The 
prevalence reported in several studies indicated that asthma may 
still be increasing in western countries and developing countries 
[2-4]. Asthma is a complex disease and has a strong genetic 
component in its pathogenesis. So far, considerable efforts have 
been made to evaluate the association between genetic variants 
and asthma risk, and numerous genes have been identified as 
asthma susceptible genes [5-7]. 

The P2-adrenoceptor (P2AR) mediates the physiological 
responses in the airway, which include bronchodilation, bronch- 
oprotection, enhanced mucociliary clearance [8] . The P2AR gene 
is located on chromosome 5q31-q32, a region that is genetically 
linked to asthma and related phenotypes [9] . There are three best- 
known polymorphisms in the coding region of the P2AR gene that 
can modulate the function of the receptor [10]. Many studies have 



investigated these polymorphisms to assess their potential contri- 
butions to the risk of asthma. It was observed in some studies that 
P2AR gene may play a significant role in the pathogenesis of 
asthma [11—16]. However, other studies showed no association 
between asthma and P2AR gene polymorphisms [17-22]. A large 
population-based study and three meta-analyses have also shown 
inconsistent results [23-26]. Since then, additional many studies 
with large sample sizes about P2AR polymorphisms on asthma risk 
have been reported. 

Therefore, we present the results of a comprehensively updated 
meta-analysis of all relevant published data to investigate the 
association between P2AR gene polymorphisms and asthma risk 
with focus on Arg/Glyl6 and Gln/Glu27 polymorphism. 

Materials and Methods 

Publication search 

We followed criteria of the Preferred Reporting Items for 
Systematic Reviews and Meta-analysis (PRISMA) to perform this 
meta-analysis (Dox SI). Two investigators (H.X. and Y.H.) 
independently carried out a systematic search of the PubMed, 
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EMBASE, BIOSIS Previews databases for articles published until 
3 January, 2014 to identify potentially relevant articles. The search 
strategy utilized in our study were as follows: "asthma or bronchial 
asthma or bronchial hyperreactivity or allergy or atopy" and 
"ADRB2 or (52-adrenergic or beta adrenergic or adrenergic or 
adrenoceptor or beta2 AR" and "polymorphism* or variant* or 
genetic* or mutant*". No publication language restrictions were 
imposed. All the searchable studies were retrieved, and we also 
checked their references for other relevant publications. The 
detailed search strategy for this study is showed in Dox S2. 

Inclusion and Exclusion Criteria 

The inclusion criteria of our study were as follows: (1) Human 
studies evaluating associations between (3 2 -Adrenoceptor gene 
polymorphisms and asthma risk, (2) There were at least two 
comparison groups, for example, asthma versus control (non- 
asthma) group, and (3) Genotype distributions in comparison 
groups should be available for estimating an odds ratio (OR) with 
95% confidence interval (CI). Studies would be discarded if they 
met one of the following criteria: (1) Not relevant to P 2 " 
Adrenoceptor polymorphisms or asthma risk, (2) The design 
based on family or sibling pairs, (3) Reviews or abstracts without 
useful information. If the same patient population was reported in 
several publications, we included the most complete study in our 
meta-analysis. If original genotype frequency data were not 
reported, we sent an email to the corresponding author for 
additional data. Studies were excluded from our meta-analysis if 
their authors did not provide us with related data. 

Data Extraction 

In the data collection process, two investigators (H.X. and Y.H.) 
independendy examined full manuscripts of eligible studies, and 
relevant data were extracted into predesigned data collection form. 
We verified accuracy of data by comparing collection forms from 
each investigator. Any disagreement was resolved by discussion, or 
a third author (Y.C.) would evaluate these articles. The following 
information was collected from each eligible study: first author's 
name, year of publication, original country, ethnicity, sample size, 
genotyping method, age, and genotype numbers in cases and 
controls. 

Publication Bias 

Publication bias of studies was assessed using funnel plot, and 
P<0.05 was considered statistically significant. In addition, 
publication bias was also evaluated via Egger's linear regression 
test. 

Statistical analysis 

We tested the Hardy- Weinberg equilibrium (HWE) in control 
group using the Chi-square test. The pooled ORs were calculated 
for allelic genetic model, additive genetic model, dominant genetic 
model and recessive genetic model, respectively. Subgroup 
analyses were performed in term of ethnicity and age. The 
heterogeneity between the studies was investigated by using 
Cochrane Q_-test. I 2 was also used to test the heterogeneity among 
the included studies. A p value>0. 10 for the Q_-test indicates a lack 
of heterogeneity among the studies, then the pooled OR estimate 
of each study was calculated by the fixed effect model. Otherwise, 
the random effect model was used. 

Statistical analyses were performed using the Revman5.2 
software (Nordic Cochrane Center, Copenhagen, Denmark), 
STATA 12.0 software (Stata Corporation, College Station, TX). 



A p value<0.05 was considered statistically significant, except for 
test of heterogeneity where a level of p value<0.10 was used. 

Results 

Characteristics of studies included in the Meta-analysis 

The flow chart in Figure 1 outiined the study selection process. 
A total of 2205 articles were identified at the initial search. 730 
duplicated studies were eliminated first. After reading the abstracts 
and tides, 1320 articles were subsequently excluded. The 
remaining 155 articles were then assessed for inclusion and 120 
articles were excluded (98 records excluded for irrelevance, 15 
records no detailed genotypes, 4 depart from HWE, 3 duplicate 
publications). In this process, study references were checked as well 
and 2 relevant publication were identified [20,27]. Finally, a total 
of 37 studies met the inclusion criteria, and were included in the 
meta-analysis [11-22,27—51]. These eligible studies contained 
6648 asthma cases and 15943 controls, and 11 Caucasian, 18 
Asian and 8 other populations. The distributions of the (32 AR gene 
polymorphism in normal controls were consistent with HWE in 
every study. Details of each study included in this meta-analysis 
were summarized in Table 1. 

(32AR Arg/Gly16 and Asthma Risk 

Among 37 studies included in the meta-analysis, 29 investigated 
the contribution of (32 AR Arg/Glyl6 polymorphism to asthma 
risk (Table SI). Overall, slighdy significant associations were found 
in allelic genetic model (OR= 1.06, 95% CI = 1.01- 1.12, 
p = 0.02), additive genetic model (OR=l.ll, 95% 
CI= 1.02-1.21, p = 0.02) for Arg/Glyl6. Stratified by ethnicity, 
significant associations were also found in Asian population in 
allelic genetic model (OR= 1.14, 95% CI = 1.05- 1.23, p = 0.001), 
recessive genetic model (OR = 1.26, 95% CI = 1.12-1.43, 
p = 008) and additive model (OR= 1.23, 95% CI = 1.06-1.44, 
p = 0001), but not in Caucasian population. No association was 
found in adults and children subgroup in term of age (Table 2). 

(32AR Gln/Glu27 and Asthma Risk 

There were 22 studies investigating the association between 
Gln/Glu27 polymorphism and asthma risk (Table S2). One study 
[32] may be the main contributor to heterogeneity and has 
significandy influenced pooled results (discuss below). Therefore, it 
was not included in data synthesis process. Overall, no significant 
association was found in all genetic models when combined all 
other 21 studies (Table 3). Age stratification showed significant 
associations in adults in allelic genetic model (OR = 1 . 1 0, 95 % 
CI = 1 .0 1 - 1 . 1 9, p = 0.03) and additive genetic model (OR =1.20, 
95% CI = 1 .07 — 1 .34, p = 0.002), but no association was found in 
ethnicity stratification (Table 3). 

Heterogeneity exploration and sensitivity analysis 

For (32 AR Arg/Glyl6 polymorphism, there were statistically 
significant heterogeneities in all genetic models when all eligible 
studies were combined (Table 2). We thus performed subgroup 
analysis on the basis of age and ethnicity. Statistically significant 
heterogeneity remained shown in Asian and Caucasian population 
subgroup, children subgroup and adults subgroup (Table 2). 
Therefore,we used Galbraith plots to graphically evaluate the 
source of heterogeneity. As showed in Figure 2A, four studies 
[31,35,47,50] may be the main contributor to the heterogeneity in 
Asian population subgroup, and Figure 2B indicated that one 
study [15] maybe the main contributors to the heterogeneity in 
Caucasian population subgroup. In order to assess the reliablity of 
our meta-analysis, we further performed a sensitivity analysis 
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3 duplicate publications 



37 articles included in this meta-analysis 



Figure 1 . Flow Diagram for Inclusion of Studies in Meta-analysis. The initial search identified 2205 articles, of which, 37 were included in the 
final analysis. 

doi:1 0.1 371 /journal.pone.01 01 861 .g001 



through sequentially excluded individual study and similar result 
was obtained. 

For (32AR Gln/Glu27 polymorphism, statistically significant 
heterogeneity were found in all genetic models, but quantitatively 
moderate heterogeneity with I 2 values<50% (Table 3). Sensitivity 
analysis through sequentially excluded individual study shown that 
one study [32] maybe the main contributor to heterogeneity and 
has significantly influenced pooled results. Combined results 
excluding this study showed small heterogeneity for both Asian 
and Caucasian population (Table 3). 

Publication bias 

For (32 AR Arg/Glyl6 polymorphism, as shown in Figure 3 A, 
the shape of funnel plot revealed asymmetry to some extent. For 
|32AR Gln/Glu27 polymorphism, the shape of funnel plot did not 
indicate any evidence of significant asymmetry (Figure 3B). 
However, the Egger's test yielded no evidence of publication bias 
for (32AR Arg/Glyl6 polymorphism (A vs G, p = 0.357), and Gin/ 
Glu27 polymorphism (A vs G, p = 0.764). 

Discussion 

Asthma is a public health problem, affecting approximately 300 
million people worldwide [52]. The mainstay of treatment for 
acute asthma includes inhaled short-acting |32-adrenergic receptor 
(|32AR) agonists, such as albuterol [5,53]. There is mounting 
evidence that polymorphisms in the p"2AR gene are associated 
with significant variability in response to short-acting |32-agonists 
(SABA). In adults, genetic variations in this receptor have also 
been linked to asthma severity [12,54]. In children, significant 
association was observed between favorable therapeutic response 
to inhaled |32-adrenergic agonists and the Arg/Arg phenotype at 
position 16 of the |32AR [55]. However, a large cohort study has 
shown that important (32AR gene polymorphisms are not main 
determinants of asthma incidence or prevalence in the British 
population [25]. Therefore, conflicting data regarding the role of 
(52AR polymorphisms in asthma susceptibility and presentation 
has been reported. 

We performed a systematic review and meta-analysis, including 
37 studies involving 6648 asthma patients and 15943 controls, to 



investigate the associations between (32 AR Arg/Glyl6 and Gin/ 
Glu27 polymorphisms and asthma. The combined results of 29 
studies showed that A allele carriers of Arg/ Gly 1 6 polymorphism 
is strongly associated with asthma risk in the total population. 
Stratification by ethnicity revealed statistically significant associa- 
tion in the Asian population. Nonetheless, overall we failed to find 
any association between Gln/Glu27 polymorphisms and asthma 
risk for 2 1 eligible studies. Subgroup analysis on the basis of age 
and ethnicity showed significant associations in adults in allelic 
genetic model and recessive genetic model for Gln/Glu27. 
However, no association was found in Asian and Caucasian 
population subgroup based on ethnicity under any genetic model. 
Sensitivity analyses through sequentially excluding individual 
study yielded similar results in both |32AR Arg/ Gly 16 and Gin/ 
Glu27 polymorphisms. In our study, we could not assess haplotype 
effects, which would have required genotype data that were 
unavailable in the published studies. 

So far, three previous meta-analyses have been conducted to 
explore the function of (32AR gene polymorphisms in asthma 
[23,24,26]. One of them shown that the Glyl6 allele of the (32- 
adrenergic receptor gene increases the risk of nocturnal asthma, 
but does not alter the risk of mild asthma [24] . The other reported 
that Gly 16 may have a protective effect for children and Glu27 
carriers may decrease the risk of asthma [23]. The third one 
showed that ADRB2 does not contribute to susceptibility of 
asthma [26]. Interestingly, our meta-analysis results are different 
from those of published meta-analysises, showing that Arg 16 may 
increased the risk of asthma in total population, especially for 
Asian population. For Gln/Glu27, significant associations could 
only be found in adults in allelic genetic model and recessive 
genetic model. The inconsistency between the results of ours and 
published articles may be partly attributed to the following factors. 
First of all, while in our meta-analysis 37 studies (29 for Arg/Glyl6 
and 21 for Gln/Glu27, respectively) involving 6648 asthma 
patients and 15493 controls were included, Migita,0 [26] 
included 8 studies for Arg/Glyl6 and 9 for Gln/Glu27, 
Thakkinstian,A and colleagues [23] only included 12 studies for 
Arg/Glyl6 and 12 for Gln/Glu27, Contopoulos-Ioannidis,D.G 
[24] included 17 studies for Arg/ Gly 16 and 16 for Gln/Glu27 
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Figure 2. Galbraith plots of the association between p2AR 
Arg16Gly polymorphism and asthma. A) Asian population 
subgroup; B) Caucasian population subgroup. Four studies^ maybe 
the main contributors to the heterogeneity Asian population subgroup; 
One studies# maybe the main contributors to the heterogeneity Asian 
population subgroup. 
doi:1 0.1 371 /journal.pone.01 01 861 .g002 

respectively. Second, some articles departing from HWE have 
been included in Contopoulos-Ioannidis D. G and Migita, O's 
meta-analysis [24,26]. Deviation from HWE is taken as an 
indication that the alleles are not segregating independently, either 
for genetic reasons or methodological reasons [56]. Therefore, this 
may be a potential source of inconsistency between our study and 
other meta-analyses. 

Genome-wide association study (GWAS) has been productive 
genotyping method to test a vast number of single nucleotide 
polymorphisms and assess their relations with complex diseases 
and phenotypes [57]. Since the first GWAS investigating 
susceptibility gene for asthma published in 2007 [58], more than 
ten GWAS of asthma have been performed in Caucasian, 
Mexican, and African-ancestry populations [51,59,60] and more 
recently in Japan and Singapore of Asian populations [61-63]. 
Nevertheless, neither an individual GWAS study nor meta-analysis 
of GWAS has implicated ADRB2 the susceptibility gene for 
asthma. In our meta-analysis, we found that neither ADRB2 Arg/ 
Glyl6 nor ADRB2 Gln/Glu27 has significant association with 
asthma in Caucasian population, which further confirmed the 
results of previous published single GWAS and meta-analysis 
[51,59,60]. However, significant associations were found between 
ADRB2 Arg/Glyl6 and asthma in Asian population, and this 
result contrasts with the results of studies recently conducted in 
Japan and Singapore [61-63], in which no single nucleotide 
polymorphisms in ADRB2 gene was found statistically significant. 
To date, GWAS studies for asthma have not been conducted in 
Chinese subjects, we are not sure whether the genomic findings of 
recent studies are applicable to Chinese. Therefore, in order to 
further clarify the role of ADRB2 gene in the pathogenesis of 
asthma, in the future, more GWAS should be conducted in Asian 
population, and especially Chinese considering its enormous 
population. 

Several strengths of this meta-analysis could be listed as follows. 
First, we followed a rigorous protocol of systematic review, 
identifying data from 3 different databases comprehensively. 
Second, we have conducted subgroup analysis, sensitivity analysis 
and Galbraith plots to explore the source of heterogeneity, 
suggesting the reliability of our study. However, the results of our 
meta-analysis should be interpreted with caution due to some 
unavoidable limitations. First, potential heterogeneity and con- 
founding factors may have affected the analysis. Statistically 
significant and quantitatively moderate-to-high heterogeneities 



Figure 3. Funnel plot of publication bias for P2AR polymor- 
phism under recessive genetic model. A) Arg/Gly16; B) Gln/Glu27. 
Funnel plot with pseudo 95% confidence limits was used. 
doi:1 0.1 371/journal.pone.01 01 861 .g003 

were found in all genetic model for Arg/Glyl6, even if we have 
found out the potential heterogeneity contributors. Second, 
publication bias might result in the loss of studies, which 
subsequently could affect the meta-analysis results. The shapes 
of funnel plots revealed slight asymmetry for Arg/Glyl6 and we 
should pay attention to this point. Third, asthma is a heteroge- 
neous group of conditions that result in recurrent, reversible 
bronchial obstruction [5]. As shown in Table 1 that the asthma 
definition varied among different articles and this may be a 
confounding factor. Case-control studies with clear diagnosis 
criteria and rigorous quality control are needed to conducted in 
the future. Fourth, most of the studies included in our meta- 
analysis were performed in Asians and Caucasians. Thus, our 
results may be applicable only to these ethnic groups. Fifth, 
although we have tried our best to contact the corresponding 
authors of published original papers and meta-analysis, we still 
cannot gain genotypes data of several studies. Sixth, so many tests 
have been conducted in our meta-analysis that it is not sure 
whether we have found something significant by chance in this 
process. Finally, asthma itself is a complex and multifactorial 
disease and potential interactions among gene-gene and gene- 
environment should be considered. 

Conclusion 

The present meta-analysis suggests that (32 AR Arg/Glyl6 
polymorphism may be an important genetic factor in the overall 
risk for developing asthma, especially in Asian population, and 
Gln/Glu27 polymorphisms may be a contributor to asthma 
susceptibility for adults. GWAS or more case-control studies with 
explicit diagnosis criteria and rigorous quality control are needed 
to support this finding in various ethnic groups. 

Supporting Information 

Table SI Distribution of Arg/Glji\6 genotypes among patients 
with asthma and controls included in the meta-analysis. 
(DOC) 

Table S2 Distribution of Gln/Gln27genotypes among patients 
with asthma and controls included in the meta-analysis. 
(DOC) 

Dox SI PRISMA 2009 Statement checklist. 
(DOC) 
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